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We study both analytically and numerically a coupled system of spherically symmetric SU(2) 
Yang-Mills-dilaton equation in 3+1 Minkowski space-time. It has been found that the system 
admits a hidden scale invariance which becomes transparent if a special ansatz for the dilaton field 
is used. This choice corresponds to transition to a frame rotated in the lnr — t plane at a definite 
angle. We find an infinite countable family of self-similar solutions which can be parametrized by 
the N - the number of zeros of the relevant Yang-Mills function. According to the performed linear 
perturbation analysis, the lowest solution with N = only occurred to be stable. The Cauchy 
problem has been solved numerically for a wide range of smooth finite energy initial data. It has 
been found that if the initial data exceed some threshold, the resulting solutions in a compact region 
shrinking to the origin, attain the lowest N = stable self-similar profile, which can pretend to be a 
global stable attractor in the Cauchy problem. The solutions live a finite time in a self-similar regime 
and then the unbounded growth of the second derivative of the YM function at the origin indicates 
a singularity formation, which is in agreement with the general expectations for the supercritical 
systems. 
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I. INTRODUCTION 

Singularity formation in nonlinear evolution equations 
is a general problem which arises in various branches 
of mathematical physics. The term singularity means 
that the solutions to the evolutional differential equations 
cease to be differentiated at some point/region. However, 
the strong-field regime of the approach to the singularity 
can provide some universality and additional symmetries 
in the solutions behavior that is of particular interest 
in the corresponding theoretical field models, including 
gravity. 

It is well-known the space-time singularities are the 
most generic features of Einstein's equations and it is 
believed that the mixmaster type singularity can pretend 
to be a generic one. However, the nature of the space- 
time singularity is model-dependent and various matter 
fields provide various scenario of the singularity forma- 
tion. In the typical case of a gravitational collapse the 
late-time dynamics of the space-time singularity develop- 
ment is hidden under the formed event horizon. However 
the last decade studies of the massless fields collapse ini- 
tiated in a pioneer work by M. Choptuik [2( opened a new 
direction for understanding the singularity formation dy- 
namics. 

It was realized Q that black holes with an arbitrary 
small masses would be obtained and the mass scale low 
was discovered. This was called Type-II behavior which 
is characterized by a mass gap absence in the black hole 
spectrum. The limiting case of a vanishing black hole 
mass is of particular interest. It has been found these 
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critical solutions observed at the threshold of the black 
hole formation are discretely self-similar. 

Later on, the numerical studies of the gravitational 
collapse of the self-interacting massless fields were per- 
formed 0, Q| and the Type-I behavior was observed as 
well. The Type-I behavior is characterized by a finite 
mass gap in the black hole spectrum. This Type-I behav- 
ior takes place if the considered system of Einstein-matter 
equations admits static finite energy asymptotically flat 
solutions. In this case the smallest black hole has a finite 
mass which is equal to the mass of the lowest static solu- 
tion. In the Einstein- Yang-Mills (EYM) system of equa- 
tions the N = 1 Bartnick-McKinnon (BK) static solution 
P| is the lowest mass static solution which determines a 
minimal mass of the formed black hole. Moreover, this 
N = 1 BK solution is occurred to be an intermediate 
attractor which the collapsing solution should attain in 
order to turn to the Type-I black hole formation scenario 

a, @. 

On the other hand the analysis of the blowup in the 
nonlinear wave equations in various field models without 
gravity showed the singularity formation in gravity and 
blowup in the nonlinear wave equations share many com- 
mon features Q, @, @, E3- Based on these observa- 
tions, P. Bizon and Z. Tabor put forward the conjecture 
13 that all basic properties of the gravitational collapse 
of massless fields such as universality, self-similarity and 
mass scaling, originally observed for Einstein's equations 
are just the basic properties of a wide class of a supercrit- 
ical evolution PDEs. This class includes Einstein equa- 
tions, Yang-Mills equation in 5 + 1 Minkovski space and 
many others. 

Following this conjecture, we consider a coupled sys- 
tem of Yang-Mills-dilaton (YM-dilaton) equations in 3+1 
Minkovski space which is of interest for several reasons. 
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First of all, this system is a truncated version of a theo- 
retical field model inspired by the heterotic string. Then 
the dilaton field itself which is also called a scalar gravi- 
ton provides many key features characteristic of a gravity. 
For example the static system of spherically symmetric 
SU(2) Yang-Mills-dilaton equations has a countable in- 
finite set of regular finite mass solutions which are 
similar to the BK solutions in the Einstein- Yang-Mills 
system of equations. Note that a similar set of the reg- 
ular solutions exists in the Einstein- Yang-Mills-dilaton 
system [l^, 0|- This fact stresses the dominant role 
of the Yang-Mills field for the YM-dilaton, EYM and 
EYM-dilaton systems dynamics. One can expect that 
the YM-dilaton model should exhibit all main features 
relevant for the EYM and EYM-dilaton models. Hence a 
study of the singularity formation in the YM-dilaton sys- 
tem should shed new light on a singularity development 
hidden under the formed event horizon in the EYM and 
EYM-dilaton models. 

We have found that the considered system of spher- 
ically symmetric SU(2) YM-dilaton equations admits a 
hidden scale invariance. The corresponding criticality in- 
dex is equal to +1 similar to the Einstein's equations 
in 3 + 1 dimension and Yang-Mills equation in 5 + 1 
Minkovski space-time. This means the YM-dilaton is su- 
percritical and singularity should arise if the initial data 
exceeds some threshold. Note a pure Yang-Mills system 
of equations in 3 + 1 Minkovski space-time is subcriti- 
cal and the solutions should remain everywhere smooth 
during the evolution We have found the late-time 

asymptotics of the YM-dilaton solutions prior to the sin- 
gularity formation is universal and is described by the 
lowest N = self-similar profile. This N — profile for 
YM field is occurred to be similar to those found in 0, 
Q for pure YM field in 5 + 1 Minkovski space. We have 
shown the whole family of the self-similar solutions exists 
labelled by the N = 0, 1, 2, 3 . . . oo - number of the nodes 
of the relevant YM function. 

The paper is organized as follows. In the next Sec- 
tion we introduce main definitions and discuss the scale 
properties of the obtained YM-dilaton system of spheri- 
cally symmetric equations. In the third Section the self- 
similar solutions of the YM-dilaton system and their lin- 
ear stability are considered. And in the fourth Section 
the results of the numerical simulations of the evolutional 
Cauchy problem are discussed. We conclude with the 
main results and discuss briefly some open questions in 
the last Section. 



II. MAIN EQUATIONS AND DEFINITIONS 

We consider a coupled system of Yang-Mills-dilaton 
fields, which is given by the action: 



S 



1 

-in 



(9$) 2 



exp (£;<!>) 



d 3 xdt, (1) 



2 X ~~' V 
where $ — dilaton field and F afiu — Yang-Mills field. 
Note that this action is a truncated version of the bosonic 
part of the heterotic string effective action in four dimen- 



sions in Einstein frame |15| . 

In the case of the spherical symmetry, the dilaton field 
and the SU(2) Yang-Mills potential can be parameterized 
in terms of two independent functions r) and f(t,r) 
as follows: 

$ = *(r,i), — , A? = e aik ^ (/(r,i)-l). (2) 

After substitution of this ansatz into the action and 
to the field equations, the rescaling <!> — > $/fc, r — > 
(k/g)r, t — * [k/g)t and S — > g 2 S removes the depen- 
dence on k and g and therefore we can put k = 1, g = 1 
in what follows without restrictions. After integrating 
out the angular variables, the resulting action and the 
field equations become: 



I r 2 $' 2 _ I r 2$2 
2 2 



f + f$-f"-f'*' 



2$' 

r 



2r 2 



/(I-/ 2 ) 



dr dt, 



(f 2 If 
2r 2 



(3) 



(4) 



.(5) 



Here and below (except Eq. (|26|l ) dot stands for a par- 
tial derivative with respect to t, while prime is a partial 
derivative with respect to r. 

It is also useful to write down the expression for the 
total energy 



E = 



2 2 



a . i2 , (/ 2 -l) 2 



.r + / 2 + 



2r 2 



dr, 



(6) 



which is, of course, conserved on shell. 

For further analysis it is important to note that the 
system Eqs. (0J 01 admits a hidden scale-invariant form. 
Indeed, the Eq. Q is scale-invariant in sense that if 
f(t, r), r) is a solution to equation @, then f(t, r) — 
f(j, j), <!>(£, r) = j) is also a solution. The same is 
not true for Eq. © because of the factor e*/r 2 on the 
right-hand side of equation JSJ, which breaks the scale 
invariance. However, it is possible to extract a scale in- 
variant part (f>(t, r) from the dilaton function <i>(i, r) that 
makes transparent the hidden scale invariance of the sys- 
tem. Indeed, if we consider the following ansatz 



$(r,t) = <f>(r,t) + 21nr, 



(7) 



the system Eqs. (JUISJl, rewritten in terms of f(t, r), <j>(t, r) 
becomes scale-invariant. The energy, expressed in terms 
of the functions f(t,r), <fi(t,r) becomes 
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providing the corresponding homogeneous scale low for 
the energy as = AW (t, r), <f>(t, r)). Ac- 

cording to the PDE general theory [16j the degree a of the 
scale parameter A as it enters in the energy homogeneous 
scale low defines the criticality class of the PDE system. 
Criticality class of the system indicates the possibility of 
a singularity formation in the corresponding well-posed 
Cauchy problem as follows. If a < the system is sub- 
critical, then all initially regular solutions should remain 
globally regular during the evolution. If a > the system 
is supercritical and one should expect singularity forma- 
tions for a finite time for all initial data, if its exceed a 
some threshold values. If a — 0, the system is critical 
and there are no definite expectations on singularity for- 
mations. Since in our case we have a — +1, the system 
is supercritical and one should expect singularity forma- 
tions, which will be confirmed in the last section of this 
paper. 

The revealed invariance of the solutions under the scale 
dilations f(t,r) — > f(j,j),(j)(t,r) — > <fi(j,j) means that 
both, t and r variables enter effectively the solutions only 
in particular combinations - r/t, or t/r. Because of the 
time translation invariance it is useful to introduce some 
positive constant T, which transforms similarly under the 
dilations, and search for general scale-invariant solutions 
in a self-similar form: 



f(r, t) = f(x), cf)(r, t) = cf)(x), x 



(9) 



The constant T has the meaning of a blowup time 
- absolute value of the time in the evolution Cauchy 
problem, when the expected singularity starts develop- 
ment and the solution cease to be smooth. We also 
will use alternatively the inverse independent variable 
r\ = i = which is more natural for study of the 

Cauchy problem. The coordinate r\ covers a half of the 
complete Minkowski space-time only, corresponding to 
the past region of the blowup point t = T, r = 0. The co- 
ordinate x covers complete Minkowski space and we will 
use it mainly for the analysis of the self-similar solutions 
in the next section. 



As we will see below, s(x) represents the regular part 
of the function <fi(x) on the semiaxis x S [1, oo). 

In terms of the functions f(x) and s(x) the system of 
PDEs Q] EJ) transforms to a system of ODEs: 
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(10) 
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This system has 4 singular points x = — oo, —1, +1, 
+oo. At the first step we restrict ourselves to the interval 
x € [l,+oo) which covers the interior of the past light 
cone of the point t = T, r = 0. 

The natural requirement of the regularity for the YM 
function fix) on the past light cone of the point t — 
T, r = provides the following local solution of the sys- 
tem (fTUl Ull near x = 1: 

/(z)x-i = h( x ~ 1) - j (1° + s i)( x ~ I) 2 + 

+f 1 (—- -f 2 + — s \ + — Sl ) (x - If + OHx - l) 4 ), 
■ /1 V48 3' n 96 1 32 V ; u ' h 



s(x) x ->i =4 + si(x - 1) + 



(12) 
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96 Sl ^ 48 1 3 



ix-lf 



+0((x-l) 4 ), 



where f\ and s\ are free parameters. The regularity re- 
quirement of the YM function / at the origin on each 
slice t < T (until the blowup time t = T) leads to the 
following series expansion of the solutions near x = +oo, 
written by making use of the inverse r\ — 1/x self-similar 
variable: 



III. SELF-SIMILAR SOLUTIONS AND THEIR 
LINEAR STABILITY ANALYSIS 

In order to bring the system of the scale invariant equa- 
tions to the form, suitable for further analysis, it is con- 
venient to introduce a new function six) as follows: 

4>ix) = In (x 2 six)) . 



fi V ) 7 ^0 = ±l + f2V 2 



= s + soisofi - ^)v 2 + 



-so i&stfl + 8s /| -q)V* + Oirf). 



(13) 
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FIG. 1: Self-similar solutions for N = 0,1, 2, 3: YM functions 
/jv(x) - top, dilaton functions sn{x) - bottom. Note, as N 
increases the solutions sn(x) tends uniformly to the limiting 
curve -4-. 



Here f 2 and sq are also free parameters. 

These local solutions at the singular points satisfy ap- 
propriate boundary conditions: 



/(1) = 0, «(1)=4, /(oo)=±l, 
f, x (x) -> , s, x (x) -* . 



(14) 



As a result we get a suitable desingularised boundary- 
value problem (BVP) for the system (|10l with the 
boundary conditions of Eq. 114(1. Alternative strategy 
consists in solving of the system l(l()l II 111 as an initial 
value problem starting with the series expansion l|12|) 
near the x = 1. In this case the problem transforms 
to the shooting method search for some discrete sets of 
free parameters /i,si, so that they provide the meet- 
ing of the local solutions to Eq. I|13|) at infinity. The 
method of shooting technique also gives a powerful tool 
to prove global existence, developed by Breitenlohncr, 
Forgach and Maison [lj} . We are planing to consider the 
proof of the existence in a separate paper. 



In order to find the above solutions numerically, we 
consider the above BVP using the method, described in 
details in ^^|. In accordance with the general expecta- 
tions, we have found that the BVP fTTH), dJ, fH|) admits 
infinitely many solutions and all of then can be labeled 
by the N - the total number of zeros of the YM func- 
tion /(x) on the open semiaxis x € (1, oo). The solutions 
with the lowest N, N = 0, 1, 2, 3 are given in Fig.^ One 
can also extract the values of the corresponding free pa- 
rameters /i, Si, or f 2 , Sq as they enter the local series 
expansions Eq. l(P2*)l and Eq. (|T3)l . The values of fi and 
Si are displayed in Tabled 



TABLE I: The values of free parameters /i, si for the self- 
similar solutions with various N 



N 



h 



0.498934096775465 
1.1357110" 12 
8.9871810~ 13 
8.0687110" 

0.0 



-r3 



-8.92179247 
-8.00124 
-8.00095 
-8.00068 

-8.0 



Note that the solution which corresponds to N = oo 
is simply the trivial solution /(x) = 0, s(x) = 4/x 2 , that 
in terms of the primary function reads 4>(x) = In 4. 

If N is finite and N > > 1 , the solution of the system 
of Eqs. (jTUl [n\i can be easily found by applying the lin- 
ear perturbations theory around the trivial solution in 
the oscillation region. Indeed, after the substitution of 
/(x) = + e/(x), 0(x) = In 4 + e<j){x) into the system 
Eqs. (|10llll|) . the linearized system for the perturbations 
f(x),4>(x) which holds at the interval 1 <^ xo < x < x^ 
(xn ~ exp(27r/v / 3A r )), becomes 



dx 2 x 2 ' dx 2 x 2 



(15) 



The solution, bounded at infinity can be found in terms 
of elementary functions as follows: 




9i 

X 



(16) 



where C\, C 2 and 8 are some integration constants. This 
oscillating solution gives a strong argument in favor of 
the existence of solutions with arbitrary N > zeros of 
YM function. 

We found the solutions, which are defined on the in- 
terval x G [1, +oo). Now one can analytically continue 
them in terms of 

the functions /(x), </>(x) to the left of the point x = 1. 
It has been done numerically. The continued functions 
behave monotonically at the interval x £ [0, 1] and they 
meet their corresponding values at the regular point x = 
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0: f{x = 0) = foe, <fc{% = 0) = <f>oo which are also their 
asymptotic values at the spatial infinity r — > +oo on each 
slice t = const, t < T in the self-similar regime. This 
can be seen from the series expansion of the solutions in 
terms of the r\ variable at r\ — > +oo (/oo, 4>ooi fx, 01 are 
free parameters): 

f(v)r,^oo = foe + /l^T 1 + 

+ ^ {-Mi + foo- fH) n- 2 + 0(n- 3 ), 



1 - ■* 'f 2 _ 1 i f 2 _ 1 f 4 



(17) 

7 ? - 2 +0(^ 3 ). 



For example, the values /(a; = 0) = = 
-0.5072593..., <j>(x = 0) = 4>oo = 2.1214115... corre- 
spond to the N — solution (see below Fig.|SJ). 

Further continuation of the solutions to the region of 
the negative x meets an obstacle at the next singular 
point x = — 1. It was found that all self-similar so- 



lutions (except the trivial one), labeled by the N — 
0, 1, 2, 3, . . . oo ceased to be differentiable in the vicinity 
of the singular point x = — 1 and therefore they cannot 
be continued smoothly to the complete manifold, covered 
by the x coordinate. 

The stability analysis of the obtained self-similar so- 
lutions can be done in terms of the linear perturbation 
theory. The negative energy states in the spectrum of the 
linearized system are usually related to the instability of 
the background solutions. 

Now we consider the problem on the semiaxis x £ 
[1, +oo). Let us introduce in addition to the variable x a 
second independent variable r so that the set of the lines 
x = const is orthogonal to the set of the lines r = const: 



T-t 



t = - In y/(T - t) 2 - r 2 . 



(18) 



The action Q in the x, r variables after the substitution 
of the anzats Q for the dilaton field becomes: 



J 



S 




dx 



(x 2 - l) 2 \d 



(/ 2 - 1) 2 

2(x 2 - 1) 



2x 



1 \ dx 



r 



(19) 



We look for spherically symmetric perturbed solutions of the following form: 
I 



f(x,r) = f N {x) + e v / 2{x 2 - l)v{x)c luJT , <f>(x,r) = 4> N {x) + e e - * ^ 2 \J x 2 - lu(x)e lUT , (20) 

I 



where the background fpf(x),4>N(x) are some solutions 
of the BVP (Unj, EJ, G3> parametrized by the N - the 
number of zeros of the YM function f(x). It is helpful to 
use the variable p instead of the x: 



H 2 \x+l 



l + e 2p 
l-e 2 P' 



(21) 



Then, following the lines of the Ref.jllJ, one can bring 
the effective action for the perturbations to the standard 
form of two-channel Schrodingcr equation 

S=- f drdpc 2T [(x) + p (x),p + X + U X -(u 2 -l)x + x\, 



where the column \ an d t ne function A(p) are introduced 



as follows (<T2 is Pauli matrix): 



A{ P ) 



p exp(0 Ar (p)/2)/ jV (p), p 
V2 



dp. (23) 



Finally we get a matrix equation, which describes the 
spectrum of the linear perturbations 



-X,p,p + U X = n 2 X , 2 =^ 2 -l, (24) 



where 



V = \ Un Ul2 ^) = e-' lrT2A V = e l,72A 
U21 U22 



(25) 
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FIG. 2: Matrix elements Uij(p) for the background self- 
similar solution with N = 0: Un - top, U12 = U21 - middle, 
U22 - bottom 



and the elements of the matrix V are 
\ 2 sinn p 

v " = v ^ = 7d^( sMi2 '' e '"'" f '''J (26) 

2 4 sinn p 

Note that prime in formula (|26|l stands for derivative 
with respect to p variable. 

For the background solution /q (p) , 4>o (p) , which corre- 



sponds to the N = 0, the matrix elements Uu(p), Uvi{p), 
Vii{p) are shown in Fig. [5] Background solutions with 
other N > provide the matrix elements Uij of a similar 
shape. 

Using the method of the phase functions shift, in- 
troduced by F. Calogero and developed by A. De- 
gasperis [23, we have found that the self-similar solu- 
tions Jn{p)i 4>n{p) with N zeros of YM function / have 
N unstable modes in the spectrum of the linear pertur- 
bations. Hence the only self-similar solution with N = 
is linearly stable. 



IV. THE CAUCHY PROBLEM 

In this section we consider the Cauchy problem for 
the system of equations (@J EJ) starting with regular ini- 
tial data. The goal is to study the behavior of the solu- 
tions near the origin, their possible attaining to some self- 
similar solutions and their late-time asymptotics prior to 
the expected blowup. 

The system of equations (@J |SJ) is a coupled system of 
nonlinear wave equations. The boundary conditions at 
the origin and at the radial infinity as well as the initial 
conditions on the slice t = have to be defined in order 
to get a well-posed Cauchy problem. In fact, such a prob- 
lem with the initial and boundary conditions imposed is 
called a mixed type Cauchy problem. However, we will 
call it Cauchy problem, for short. 

The system of equations (0J[SJ), similarly to its static 
version, have two singular boundary points r = and 
r = +00. The regularity requirement provides the only 
possible series expansion near the origin: 

f(t,r) r ^ Q = l-b(t)r 2 +0(r 4 ), 

$(t, r) r ^ = *o(«) + <5>2{t)r 2 + 0(r 4 ). (27) 

where bit), &o(t), $2(0 are finite smooth functions. 
Hence the corresponding boundary conditions at the ori- 
gin, which must hold during the evolution, are as follows: 

/(t,r = 0) = l, /'(t,r = 0) = 0, 

$(t, r = 0) = $o(*), * '(*, r = 0) = 0. (28) 

It is important to note that the value of the dilaton 
function at the origin <£>(£, r = 0) is a free parameter and 
there are no reasons to keep it equal to its initial value. 
Indeed, in our numerical algorithm we have put only the 
condition &(t = 0, r = 0) = while $(£, r = 0) has been 
calculated on each slice t > according to the evolution 
equations I0J[SJ). This means we have put free boundary 
condition for $(i,r = 0). 

The asymptotic behavior of the solutions regular at 
infinity is given by the series expansion: 



f(t,r) 



±(l--+0(r- 2 )), 



r) r ^oo = $oc - - + 0(r~ 4 ), (29) 
r 
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FIG. 3: The set of YM profiles f(t,r) at the various times 
prior to the singularity formation. The ingoing waves at 
the small radial scales exhibit universal self-similar behav- 
ior whereas the outgoing waves at the right side of the figure 
are almost frozen at the given time scale. The solutions cor- 
respond to the following initial data data: 



I. /(0,r) = (l-ar J )/(lW)- 

II. / (0, r) = 1 - Ar 2 exp(-cr(r 
A = 0.2, a = 10, r = 2 (solid). 



kink-type, a = 0.281 (dashed); 



ro) 



Gauss- type, 



where c, d and <E>oo are constants. This gives the bound- 
ary conditions at the spatial infinity as follows: 



/(t,r = oo) = ±l, 
$(*,r = oo) = $oo 



/'(i,r = oo) = 0, 
$'(i,r = oo) = 0. (30) 



In order to impose the initial conditions in the Cauchy 
problem for the system Eqs. JU EJl, it is enough to set 
the initial distribution and its derivatives for the YM 
function f(t = 0,r) only. Then the dilaton distribution 
<3>(< = 0, r) can be obtained from the field equation J5j, 
in a way similar to the Einstein equations, where the 
dilaton plays a role of the relevant metric function now. 

So, we always consider smooth regular initial distribu- 
tions f(t = 0, r) = h(r), where h(r) provides non-vacuum 
values for the YM field function in some compact region 
for r outside the origin. It is convenient to use two dif- 
ferent types of the initial profiles for h{r): 



h(r) = 1 - A r 2 exp [-a(r - R) 2 ] 



(31) 



which is a Gauss-type (A, a and R - constants) initial 
profile that connects the same YM vacua / = +1, and 



h(r) 



I — ar 2 
1 + a r 2 



(32) 



which is a kink- type (a — const.), that connects two 
topologically distinct YM vacua / = ±1. 

After the initial distribution h(r) is fixed, one can de- 
fine YM radial wave, propagating towards the origin (in- 
going wave) as follows: 



/(0, r) = h{r), /(0, r) = f '(0, r) = h'(r). (33) 
We also put: 



$(0,r) = 0. 



(34) 



The dilaton function at t — can be obtained now from 
the Eq.JSJ) with the initial data at the origin imposed: 



2$' 

r 



1 



2r 2 



$(0,r = 0) = $'(0,r = 0) = 0. 



(35) 



The symmetric initial conditions that lead to two YM 
radial waves, propagating towards (ingoing) and out- 
wards (outgoing) the origin are determined in a similar 
way: 

f(0,r) = h(r), /(0,r) = 0, $(0,r)=0, 



2$' 



,21 



/ 2 



r 

$(0,r = 0) = $'(0,r = 0) 



2r 2 



0. 



(36) 



The Cauchy problem has been studied numerically with 
the help of a finite difference scheme that preserves the 
total energy during the evolution. Adaptive mesh refine- 
ment algorithm was also installed in order to get a better 
resolution on a small scales prior to the singularity for- 
mation. 

The system has been solved for various initial profiles 
h(r) and initial data (|35l - I36fl . The results can be sum- 
marized as follows. If we start with small initial data, 
controlled by the values of the parameters in Eqs. I|31|l . 
B32JI . the ingoing wave solution get smoothly bounced 
near the origin and after which is radiated away quali- 
tatively similarly to the typical solution of a linear wave 
equation. If we gradually change the relevant parame- 
ters in the initial conditions towards more strength initial 
data, after it exceeds some threshold value, the ingoing 
YM wave solution demonstrates universal behavior at the 
small scales. In fact, the YM function f(t, r) attaines the 
self-similar N = solution fo(jr-z), see Figs. |3 0] As 
it was expected, in this regime the corresponding solu- 
tion for the dilaton function <£>(£, r) develops in a non 
self-similar way. However, as it can be easily seen from 
Fig. the whole <&(t, r) profile propagates along the 
straight line 21nr + 2. This means that in the appropri- 
ately shifted and rotated at angle (3 = arctg(— 2) back 
clockwise frame the obtained dilaton function <p(t,r), 
<p(t,r) — $(t,r) cos /3+sin/3 In (r/r ), r = 2, exhibits the 
self-similar behavior and attains the self-similar N = 
solution <fio(t,r). The dilaton function $(t, r) can be ex- 
pressed now in terms of its self-similar part <I>(t, r) — 
4>(t,r) + 2 bar that reproduces the ansatz Q. Fig. [3] il- 
lustrates the late-time evolution of some typical solution 
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ln(r) ln(r) 

FIG. 4: Self-similar regime of the solutions behavior: YM function f(t,r) (left), dilaton function $(i,r) (right) at the various 
t prior to the blowup. The solutions effectively depends on T — t, where T is blowup time. The initial data is the same as in 
Fig. [3 kink-type initial profile (cycles) , Gauss-type initial profile (squares) . 



to the Cauchy problem which attains the self-similar so- 
lution /o(j?),^o(»?)- 

The solutions evolve in a self-similar regime during 
a finite time and at some blowup time T the second 
derivative of the YM function / exhibits unboundedly 
growth at the origin. This is an indication for a singu- 
larity formation. Blowup time T is just the total time 
in the Cauchy problem and it depends on the initial 
data. However, as the solutions attain the self-similar 
profile /o (77) , 0o (v) their dependence on the time enters 
effectively only in the form T — t. Hence, the late-time 
asymptotics becomes universal for an arbitrary solution 



with the initial data, which leads to the blowup. Accord- 
ing to our studies, the self-similar solution /o(?7), </>o(??) 
is linearly stable and, as a result, it can pretend to be a 
global stable attractor in the Cauchy problem. The sin- 
gularity formation at the origin is not accompanied by an 
energy concentration which is typical for the supercritical 
systems. Indeed, the total energy Eq. (JHJ) at time t < T 
inside the past light cone J Q t dr (. . .) of the point (T, 0) 
can be calculated using the N = self-similar solution 
fo(x), 4>o{x) as follows: 



00 , 

M = E={T-t) J dxl ^±1^{ X )1 

1 ^ 



+ 1 



fo(x)% 



(fo(x) 2 



2x 2 



4(T-t), 



and vanishes as t — > T. Note that numerically the inte- 
gral is equal to 4 with accuracy 10 -9 . The corresponding 
Schwarzschildean radius is equal to r s ~ 8(T — t). It 
means that if the gravity is included, the obtained late- 
time self-similar attractor (if exists) has to be hidden 
under the formed event horizon. 



V. CONCLUSIONS AND DISCUSSION 

We conclude with the following summary. Using the 
special ansatz for the dilaton field, we brought the system 
of the SU(2) spherically symmetric Yang-Mills-dilaton 



equations to a scale-invariant form and found an infinite 
countable family of self-similar solutions, labeled by the 
N > - number of zeros of the relevant YM function. 
Among them the only lowest solution, which corresponds 
to the N = is stable in the framework of the linear per- 
turbation theory. Being a scale- invariant, the considered 
system has a criticality index equal to +1 which means 
that the system is a supercritical one in the PDE ter- 
minology and should exhibit a singularity formation if 
the initial data in the evolutional Cauchy problem ex- 
ceeds some threshold value. The Cauchy problem has 
been solved numerically for a wide range of smooth finite 
energy initial data and the results we found are in agree- 
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FIG. 5: The late-time evolution of a some typical Cauchy problem solution prior to the blowup. YM profiles f(r/, i) (left) and 
the scale-invariant part of the dilaton profiles <j>(ri, t) = $(r, t) — 2 ln(r) — (right) attain the self-similar solution N = (dashed 
line). The values /oo,</>oo are the asymptotical values of the self-similar N = solution at r\ — > oo. The initial conditions are 
the same as in Fig. |3 for Gauss- type initial profile. 



merit with this general expectation. It has been found 
that if the initial data exceeds some threshold, the re- 
sulting solutions in a compact region, shrinking to the 
origin, attain the lowest N — stable self-similar pro- 
file, which plays the role of an attractor in the Cauchy 
problem. If the solutions attain the N = self-similar 
attractor they evolve in an universal way during a finite 
time until the second derivative of the YM function at 
the origin starts growing infinitely. 

The problem of the threshold, that separates the dis- 
perse and blowup behavior of the Cauchy problem solu- 
tions requires more detailed studies. This threshold is 
usually related to the intermediate attractors, which are 
some local minima or saddle points of the effective action. 
In our system both static N = 1 an d self-similar solu- 
tion /i (77) , (j>i (77) are saddle points of the effective action 
and can play the role of such intermediate attractors. 
The study of the decay of the static and the self-similar 
solutions with N > 1 along their unstable modes are of 
particular interest for the threshold understanding. 



As it was noted earlier, the system of Yang-Mills- 
dilaton fields is the bosonic part of the appropriately 
truncated supersymmetric field theory which also con- 
tains spin 1/2 dilatino and gaugino fields in the fermionic 
sector. The study of the fermionic fields in the obtained 
self-similar bosonic backgrounds seems to be a very in- 
teresting task itself, and also can shed new light on the 
singularity formation thresholds of the bosonic fields. 

All these tasks are under considerations and will be 
reported elsewhere. 
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